Valley splitting of AlAs two-dimensional electrons in a perpendicular magnetic field 
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By measuring the angles at which the Landau levels overlap in tilted magnetic fields (the coin- 
cidence method), we determine the splitting of the conduction-band valleys in high-mobility two- 
dimensional (2D) electrons confined to AlAs quantum wells. The data reveal that, while the valleys 
are nearly degenerate in the absence of magnetic field, they split as a function of perpendicular 
magnetic field. The splitting appears to depend primarily on the magnitude of the perpendicular 
component of the magnetic field, suggesting electron-electron interaction as its origin. 

PACS numbers: 73.21, 73.23, 73.43 



In many semiconductors, the energy dispersion of the 
conduction band contains more than one minimum, or 
valley. Examples include Si and AlAs, where electrons 
occupy pockets near or at the equivalent X-points of the 
Brillouin zone. One of the long-standing and controver- 
sial problems in the physics of two-dimensional electron 
systems (2DESs) in such semiconductors has been the 
nature of valley splitting, i.e., the lifting of this valley 
degeneracy. While it is clear that mechanisms that break 
the symmetry of the crystal potential, such as uniaxial 
strain, can lift the degeneracy, it has remained controver- 
sial whether interaction between the electrons can also 
lead to a splitting [Q. In this Letter, we report measure- 
ments of transport in low disorder 2DESs in modulation- 
doped AlAs quantum wells, which provide clear evidence 
for the dependence of valley splitting on the applied per- 
pendicular magnetic field {B±). The data reveal that 
at zero magnetic field the electrons occupy two nearly- 
degenerate in-plane conduction-band valleys. With the 
application of B±, we observe a splitting of the valley 
energies that increases monotonically with B± and is es- 
sentially independent of the parallel component of the 
magnetic field (B\\). This suggests that electron-electron 
interaction is responsible for the splitting. 

We performed experiments on 2DESs confined to mod- 
ulation doped AlAs quantum wells grown by molecular 
beam epitaxy on (100) GaAs substrates. In these sam- 
ples, the AlAs quantum well is flanked by undoped and 
Si-doped layers of Al0.4Ga0.eAs [||. We studied three 
samples with two different quantum well widths: 150 A 
for SI and S3, and HOA for S2. In bulk AlAs, electrons 
occupy three valleys that are located at X-point (rather 
than at A-point as in Si) of the Brillouin zone. Previ- 
ous studies H H H7 U have indicated that strain due to 
lattice mismatch between AlAs quantum well and GaAs 
substrate leads to a lifting in energy and depopulation of 
the out-of-plane valley in AlAs quantum wells wider than 
60A. Consistent with these studies, the 2DES in our sam- 
ples occupies two valleys whose principal axes lie in the 
2D plane. In these remaining valleys, the cyclotron elec- 
tron effective mass m* is 0.46 (in units of the free electron 
mass, mo) and the band g-factor <?{, is 2. We measured 



the magnetoresistance of the samples in L-shaped Hall 
bars aligned with [001] and [010] directions. Using illu- 
mination and front/back gate biasing, we were able to 
vary the 2D electron density n between 4.0 x 10 11 and 
9.6 x 10 11 cm~ 2 . At 30mK, the electron mobility of our 
samples was as high as 250,000 cm 2 /Vs. 

To measure the valley splitting (AEy) as a function 
of magnetic field, we utilize the coincidence method . 
This technique makes use of the difference in the mag- 
netic field dependencies of the orbital and spin energies, 
and has been widely used to measure the spin-splitting 
and the effective g-factor in 2D systems. In a 2D sys- 
tem in a strong magnetic field, quantization of the or- 
bital motion leads to the formation of Landau levels 
(LLs). The energy separation between the LLs is equal 
to the cyclotron energy, which, in an ideal 2D system 
equals heB±/(m*mo) and therefore depends only on B±. 
Thanks to the Zeeman coupling, each LL splits into two 
energy levels, one for each polarization of spin. The en- 
ergy separation between these levels, the Zeeman energy, 
is generally a function of the total magnetic field (B to t) 
and is equal to |<7*| /is B tot where /is is Bohr magneton 
and g* is the effective g-factor. In a typical coincidence 
measurement, the 2D sample is placed in a magnetic field 
whose direction makes an angle 9 with the normal to the 
2D plane. The Zeeman energy is then measured (in units 
of the cyclotron energy) from the values of 9 at which ev- 
ery other magnetoresistance minima disappear; the dis- 
appearance is the result of the overlap, at the Fermi en- 
ergy, between two energy levels associated with opposite 
spin. Now, in a two-valley system, yet another splitting 
appears: each spin-resolved level splits in two levels, sep- 
arated in energy by AEy, corresponding to each of the 
valleys. In our measurements presented here, by care- 
fully monitoring 9 at which magnetoresistance minima 
disappear, we are able to determine both the Zeeman- 
and valley-splittings. 

We first present data to establish that, at low magnetic 
fields, the 2D electrons in our samples occupy two nearly- 
degenerate valleys. Magnetoresistivity p xx as a function 
of B± (at zero B\\) for sample SI at a density of n = 9.2 x 
10 11 cm" 2 is shown in Fig. 1 |)J. Strong minima in p xx 
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FIG. 1: Magnetoresistivity p xx as a function of perpendicular 
magnetic field B± for AlAs 2D electrons in sample SI. At 
low B± , p xx in the inset shows strong minima at every fourth 
filling factor v, indicating a fourfold degeneracy of the Landau 
levels (2 for spin and 2 for valley). At higher B±, all minima 
become stronger than their (u + 4) counterpart, implying that 
both spin- and valley-splitting increase with B±. 



occur at every fourth Landau level filling factor v, which 
indicates a fourfold degeneracy of the Landau levels (two 
for spin and two for valley) JTc| ]. This degeneracy can 
be clearly seen in the Fourier transform of p xx vs. 1/B± 
data, shown in the inset to Fig. 2 for sample S2 at a 
density of 7.9 x 10 11 cm -2 . The frequency of a peak in 
the Fourier transform can be converted to a 2D density 
by multiplying by e/h, where e is the electron charge 
and h is the Planck constant. The presence of a peak 
(square symbol) at the frequency equal to one quarter of 
the frequency associated with the total density confirms 
the spin- and valley-degeneracy at low B± . We assign the 
other peaks in the Fourier transform to the spin- and/or 
valley-resolved oscillations which are observed at higher 
B± Q. 

In Fig. 2 we show a plot of the measured peak posi- 
tions in Fourier transforms as we change the density in 
S2 via a front-gate bias. The valley/spin degeneracies 
persist at all n and, as the density is decreased, all the 
peak positions linearly decrease and extrapolate to ap- 
proximately zero in the n = limit. It is particularly 
noteworthy that the low-frequency peak does not show 
any splitting within the resolution of the Fourier trans- 
form at any density; this observation indicates that the 
densities of states and therefore the effective masses of 
electrons in the two occupied valleys are the same, con- 
sistent with the 2DES occupying two nearly degenerate 
valleys whose principal axes lie in the 2D plane. 

Figure 1 reveals that, at high B±, all the p xx minima 
become stronger than their (y + 4) counterparts. For ex- 
ample, the p xx minimum at v — 9 is deeper than the 
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FIG. 2: Inset: An example of the Fourier transform of p xx vs. 
1/B± for sample S2 at n = 7.9 x 10 11 cm~ 2 . The presence 
of a peak at the frequency equal to one quarter of frequency 
associated with the total density confirms the degeneracy of 
the spin and valleys at low B± . Main: Fourier transform peak 
frequencies as a function of the front-gate bias. 



minimum at v = 13, while at v — 17 there is no visible 
minimum. This progression suggests that all the three 
relevant energies in the system, i.e., the cyclotron en- 
ergy, Zeeman energy, and AEy , increase with increasing 
B±. Our coincidence measurements, summarized in Fig. 
3(a), provide for a quantitive determination of the en- 
ergies. This figure, which is the highlight of our study, 
shows a (color) plot of p xx vs. l/cos9 (x-axis) and filling 
factor, v (y-axis). The plot was made by taking mag- 
netoresistance traces at 32 different 6>'s. A striking al- 
ternating diamond pattern emerges in the data. This 
pattern also has periodicity of four in v. Surprisingly, 
the pattern changes very little over the (B, 9) parameter 
space; as shown later, this observation implies that AEy 
increases approximately linearly with B_\_. Such a de- 
pendence leads to the simple energy fan-diagram shown 
in Fig. 3(b). Another prominent feature of the pattern 
is its right-left symmetry, which indicates that AEy is 
independent of B\ i . 

In order to make a quantitative analysis of the data, we 
assume a general model for the energies of a two- valley 
system tilted in a magnetic field: 
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where E\ and Ei are the energy levels of valleys 1 and 
2, respectively, N and N' are the LL indices, and g* 
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FIG. 3: a) Magnetoresistivity p xx vs. l/cos9 and v measured 
in sample S3 at T=30mK. Blue and red colors represent small 
and large values of p xx , respectively. To increase the visibil- 
ity of p xx oscillations over the whole field range, we scaled 
p xx by a Dingle-like factor, exp(l/6AB±). A striking alter- 
nating diamond pattern emerges in the data. This pattern 
has an approximate right-left symmetry, has a periodicity of 
four in v, and changes little over the whole field range, b) 
Fan diagram that qualitatively accounts for the data. In this 
diagram, AEy scales linearly with B± and is independent of 
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is the effective g- factor [0|. By fitting the coincidences 
between energy levels of the same valley using this model, 
we can determine the Zeeman energy (in units of the 
cyclotron energy), or equivalently the product |g*|m*, as 
a function of B± and flu. Likewise, from the coincidences 
between the energy levels of the different valleys, we can 
determine AEy (in units of the cyclotron energy), or 
equivalently, the product m* AEy. 

We first determine the product |g*|m* directly from 
coincidence angles 01,84, and 8^ in Fig. 3(a); these angles 
correspond to the crossings between the LLs associated 
with the same valley. We find that, at low fields, \g* \m* is 
constant, identical for both valleys, and at n = 9.55 x 10 11 
cm~ 2 equal to 1.76 ± 0.02. Surprisingly, while |<?*|m* is 



enhanced with respect to its band value of 0.92, it does 
not oscillate as function of v contrary to the theoretical 
expectation ||. Papadakis et al. Q have also observed 
a constant enhancement of |g*|m* in low density (n — 
2 x 10 11 cm" 2 ) AlAs 2D electrons for v > 9. In the 
higher density samples that we have studied, however, 
we find a monotonic increase in |g*|m* for v < 7. Here 
we limit our determination of AEy to the field range 
where the measured |g*|m* is constant. 

Next, we determine AEy from the values of angles 
#2, Q3, Ob and 6> 6 at which the energy levels of different val- 
leys overlap. Since we do not know exactly the absolute 
LL index, but only the difference between the indices of 
LLs to which these energy levels belong, AE V can be de- 
termined only to within an additive term 2lheBj_/m*mo, 
where I is an integer. In sample S2, the presence of coin- 
cidences at 02 for low v (not shown) requires I to be zero. 
For sample S3, coincidences indicate that the system is 
not valley-polarized for v > 8. This observation restricts 
I to either zero or one. We chose zero, since I = 1 implies 
the unlikely situation that valley splitting first shrinks to 
zero and then increases as a function of B±. 

Valley splitting measured for samples S2 and S3 is sum- 
marized in Fig. 4. Since AEy does not vary with the 
coincidence angle, it is independent of B\\ [jl3| . Further- 
more, AEy for different samples and densities fall on the 
same curve which exhibits a nearly linear dependence on 
B± for B± > 1.5T. A least squares fit of the data in this 
high-field range gives AEy = — 0.22 + 0.25J5^ (assuming 
m* = 0.46), with AEy in units of meV and B± in units 
of Tesla. At low fields, AEy deviates from the linear 
behavior. 

Two features of the data in Fig. 4 are noteworthy. 
First, since the widths of the LLs do not affect the po- 
sitions of the coincidences, there should be no disorder 
corrections to the plotted values of AEy Q. Second, 
the fact that AEy appears to be independent of density 
and and depends only on B±, suggests that electron- 
electron interaction is responsible for its enhancement 
with B±. The linear dependence on B±, on the other 
hand, is puzzling; one would normally expect a B 1 / 2 de- 
pendence as the enhancement should inversely scale with 
the magnetic length. We hope that the results presented 
in this Letter will serve as incentive for developing a the- 
ory to explain the enhancement of valley-splitting in a 
magnetic field. 
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FIG. 4: Products m*AEv and AEy (assuming m*=0.46) as 
a function of B± with typical error bars. Solid line is the least 
squares linear fit for the high field valley splitting. 
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